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ABSTRACT 



Allocation of resources has become a classical 
problem In optimization by mathematical programming. 

In the field of military applications attack aircraft 
assignment has been treated widely by deterministic 
and/or linear models. However, destruction of a target 
is no certainty nor is damage inflicted on a target 
linear with respect to the number of weapons delivered 
on it. Recent extensions J.n the field of nonlinear 
programming in conjunction with the widespread use of 
electronic digital computers permit a more realistic 
approach to this problem. This paper formulates a 
stochastic nonlinear model for assigning a force of 
attack aircraft on a single sortie against fixed location 
targets. The number of aircraft alive at weapon release 
on any pass of a series against a given target is treated 
as a random variable. The total value of damage to all 
targets is taken as the measure of effectiveness and a 
particular form of the objective function derived. The 
parameters of the model and the form of the constraint 
equations are also discussed. 
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CHAPTER I 



THE PROBLEM 
I. BACKGROUND 

Allocation of resources has become a classical 
problem in optimization by mathematical programming. 
Linear programming has been used extensively to provide 
a solution. It is valid when the relationship between 
variables is linear, however, linearity is not always 
the case. For example, an individual may find the use- 
fulness derived from owning five automobiles not nec- 
essarily five times that derived from owning Just one. 

In the past, for cases such as this, it was customary 
for one to either be satisfied with a linear approxi- 
mation or to search by other means for a solution. The 
extension of nonlinear programming over the past five 
years to its present state provides the analyst with a 
mathematical programming technique for treating the 
allocation problem with certain types of nonlinearity 
in a direct rather than approximate fashion. Dorn [l8j , 
Graves and Wolfe £l?j , Hadley ^25j , and Wolfe £l3j have 
published comprehensive works on the status of nonlinear 
programming . 
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A particular allocation problem in the field of 
military applications that has received a great deal of 
attention is that of attack aircraft assignment. Numerous 
studies and papers have been published on the subject, 
and several which have been found to be of particular 
Interest are! Limited Way Operations (U) Jioj , Tactical 
Air Warfare Force Allocat^on (U) [ 21 ] , Dresher [ 4 ], and 
Haering M. 

II. ASSUMPTIONS 

It was through participation in the Center for Naval 
Analyses Study l Tactical Air Warfare . 1964 . that the 
necessity for a formulation such as developed in this 
paper became apparent. The allocation of attack aircraft 
to fixed location targets is treated under the following 
broad assumptions J 

(Al.l)l An optimal allocation is one in' which the 
total value of damage to all targets is a maximum, subject 
to the constraints Imposed by the availability of aircraft, 
fuel, and weapons. 

(A1.2) The damage inflicted on a target is nonlinear 
with respect to the number of aircraft passes made 
against it. 

1(A identifies assumptions^ 

(D . ) identifies definitions. 
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(A1.3) There Is a probabilistic interpretation 
of the number of aircraft alive on a particular pass 
against a given target. That is, given ten aircraft 
assigned to attack a target, there is a probability 
that all ten aircraft will not be alive on any parti- 
cular pass against that target. 

III. THE FORMULATION 

Under these assumptions the allocation of attack 
aircraft was formulated as a stochastic nonlinear 
mathematical program with its attendant objective function 
and constraint equations. The objective function in this 
case is a nonlinear function of the number of attack air- 
craft alive on each pass against each target. Its 
particular form is derived and a discussion of the con- 
straint equations given in Chapter II. 

The number of passes against a target is treated as 
a random variable in the following way. the number of 
aircraft alive at release on a particular pass against 
a target is a random variable, therefore the number of 
passes against a target is a random variable which is 
the sum of the number of aircraft alive at release on 
each pass. Their probability distributions are derived 
in Chapter III and Appendix A. 

Chapter IV is devoted to a discussion of the parameters 
of the model. These include target value, fuel required, 
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racks available for ordnance, the probabilities of 
detection, engagement, kill, acquisition, and pass 
survival. 

Conclusions and recommendations comprise the final 
chapter of this paper. 
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CHAPTER II 



THE OBJECTIVE FUNCTION AND CONSTRAINT EQUATIONS 
I. THE OBJECTIVE FUNCTION 

The measure of effectiveness chosen for the 
allocation of attack aircraft was the total value of 
damage inflicted upon enemy targets. The problem is to 
maximize the total expected value of damage subject to 
certain constraints. 

(A2.1) Assume that each aircraft assigned the 
target delivers the same number of preferred weapons on 
each pass against that target. A preferred weapon 
implies that for a specific target there exists a weapon 
which is most effective in destroying that target. 

(A2.2) On any pass against a target assume that the 
target is either killed or not killed. This implies that 
a target can be destroyed on one pass. 

(A2.3) Aircraft make passes until all weapons are 
expended or until the aircraft is killed. 

(D2.1) The base-aircraft index specifies an aircraft 
type located at a specific base. For example,, 

J = 1 denotes the set of A4 aircraft from base 
number one, 

j = 2 denotes the set of F4 aircraft from base 
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number one, 



J = J - 1 denotes the set of A4 aircraft from 
base number S, 

J = J denotes the set of F4 aircraft from base 



number S. 

A. y_ 

(D2.2) denotes the number of J 1 ' base-aircraft 

a. r_ 

type assigned the i^ target, i = 1,2,..., I and 
J 1,2,... ,J. 

(D2.3) denotes the vector (N^, N i2 > ••• » 

a. y. 

(D2.4) Nj denotes the total number of the J base- 
aircraft type available. 

(D2.5) denotes the number of passes per aircraft 

planned by the base-aircraft type assigned the i 
target. 

(D2.6) PT^» a constant, denotes the probability 
the i 0 target is killed by exactly one pass given that 
b preferred weapons are delivered on that pass. 

(D2.7) (N^) denotes the random variable which 

a. y_ 

represents the total number of passes against the l 
target. W^N^) will be written as when convenient. 
The range of possible values is! 



W i - 0,1,2, ...E^N iJ R 1J . 
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(D2.8) Let be a random variable such that 
Y^ = 0 if the i*'* 1 target is not killed 

jl if the i*'* 1 target is killed at least once. 

Then 

P l * 1 = all w f ^ ’ ° I = ’] ? [“i = "] 

= E ( 1 - PT 
all w 

and 

F K = l] - 1 - p K ■ °] 

= 1 - E ( 1 - PT b ) W P f W. = wl . 
all w 1D L i J 

(D2.9) Let v^ denote the pre-assigned value of the 
th 1 

i*' targets. and let V - E v. be the total value of the 

1=1 1 




target complex. 

(D2.10) Let denote the random variable which 
represents the value of damage to the i target such 



that; 




Then, 



0 if Y t 



v i lf Y i 



0 

1. 



E [dJ = £ d‘P f D, = dl 

L 1J all d L 1 J 

-V. I 1 !- 1 ] 

■ v i i 1 - all V 1 - )W p [ M i - w ] } 

The total value of expected damage (TD) to the target 
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complex Is 



(D2 .11) TD = Z E [dJ 

11 w r 

= Z v - Z v Z ( 1 - PT ) w P [W 

1=1 1 1=1 1 all w 10 L x 

= V - Z v. £ ( 1 - PT., ) w P [w, = w 

1=1 1 all w 10 L 1 

J 

w = o^i) • • t ^ r n ^ jR ^ ^ • 

To maximize TD one must minimize the expression 

Z v 1 Z ( 1 - PT lb ) w P [w jL = w] , which 
i w 

t* Vi 

is simply the sum of the product of the value of the i^ 
target and the probability the i^ target is not killed* 
Therefore the objective function is 



(2.1) min Z v 



1 S ( 1 - n lb > 

1=1 1 all w lD 



w 



[ W 1 = B ] 



II. THE CONSTRAINT EQUATIONS 

,T f 

The minimization of the objective function is subject 
to certain constraints on the resources available. The 
four constraints given below are a non-exhaust lve subset 
of those which could be imposed, but are representative 
of the most important ones! 

1. The number of aircraft allocated to targets 
must be less than or equal to the total number of aircraft 
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available, i.e., 

(2.2) , N u * 0., 

L. V_ 

where Nj = the total number of j 6 base-aircraft type 

^ In 

available, and = the number of j base -aircraft 
type assigned the 1 target, 

2. The total fuel required must be less than 
or equal to the total fuel available, l.e., 



(2 - 3) ji f u N u- F d ■ 

where Fj = the total fuel available for the J base- 

^ Vk 

aircraft type, and f^j = the fuel required for the J 

4 “ Vi 

base-aircraft type to strike the 1^ target. 

3. The total ordnance loading must be less than 
or equal to the ordnance available, l.e., 



(2.4) 



I ' 

Z 

1=1 



b l,)n N ij 




* 



where Bj n = the number of weapons of type n available for 
the j element, and b^j n = the number of preferred 
weapons of type n the j element carries to the i 1 ' 
target. % 

4. There exists an upper limit on the number 
of aircraft which one Is willing to lose on a given 
mission, l.e., 

(2 - 5) J x L ij H id - L d • 
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+* Y\ 

where Lj = the maximum acceptable number of J base- 
aircraft losses, and = the expected percentage 
attrition of the J tJl base-aircraft type assigned the 
i th targe to 

There are numerous ramifications to these constraints 
which will not be covered in detail here c For example 
there are additional restrictions to (2.3) , (2.4), and 
to (2.5) in that, respectively, f^j is constrained by 
the amount of fuel an ( aircraft can carry, b^j n is con- 
strained by the maximum ordnance load and rack restrict- 
ions of the aircraft type, and is a function of the 

aircraft type, speed, penetration altitude, the enemy 
defenses, and the number of passes made against the 
target. 

The allocation problem has now been formulated as" 

('2.1) min Z v„ Z ( 1 - PT, K ) W P f W, = wl , 
i=l 1 all w 10 L 1 J 



subject to l 



( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 



I 

Z 

i=l 

I 

X 

i=X 

I 

z 

i=l 

I 

z 

i=l 



N iJ * N J • N iJ 



f M £ F 



b N ^ B 
D i Jn i J D Jn » 



L iJ N iJ ~ L J 



0 
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III. NONLINEAR PROGRAMS 



Rosen [ 7 ] , [ll] and Fi- acco and McCormick |^2oJ , 

[ 22 J , [^o] have developed algorithms which have been 
successful in solving nonlinear programming problems 
subject to linear or nonlinear constraints. 

The Sequential Unconstrained Minimization Technique 
(SUMT) of Flacco and McCormick and the Gradient 
Projection (GP) technique of Rosen are both available 
through the IBM Share General Program Library as Share 
Distribution 3189, RAC SUMT, and Share Distribution 1399, 
SD GP 90 respectively. 

The use of SUMT is precluded in this allocation 
problem due to the non-dif ferentiable nature of the 
objective function. That is, is a discrete integer- 
valued random variable. The concluding remarks of the 
Share write up of RAC SUMT indicate that special sub- 
routines are being developed to handle non-dif ferentiable 
functions . 

The program GP 90 should handle the non-dif ferentiable 
objective function since it uses two-sided differences 
in place of the gradients. Thus it is unnecessary to 
explicitly evaluate the gradients of the objective 
function. 
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CHAPTER III 



THE PROBABILITY DISTRIBUTIONS 
I. ASSUMPTIONS 

The probability distributions of two of the random 
variables, the number of aircraft alive at release on 
a particular pass and the total number of passes against 
a target, are presented in this chapter. 

The following assumptions were made* 

(A3.1) An individual aircraft is assigned only one 
target per sortie. 

(A3. 2) A raid is composed of one base-aircraft type, 
but more than one raid can be assigned to a target. 

(A3. 3) Enemy fighters, if scrambled against a raid, 
are sent in numbers sufficient to engage each aircraft 
in that raid. 

II. DEFINITIONS 

The following definitions were useful; 

(D3.1) PDj^n) denotes the probability that a raid 

j. Vi 

of size n is detected enroute from the j base-aircraft 

"h Vi 

location to the i target. 

(D3.2) PE^j(N,n) denotes the probability that a raid 

ty Vl 

of size n is engaged enroute from the j ' base-aircraft 

j. Vi 

location to the 1 target given that the raid is detected 
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and given that a total of N aircraft are employed in 
strike operations, 

(D3»3) PK denotes the probability that any aircraft 
in a raid is killed enroute given that the raid is 
engaged. 

(D3.4) PA^jfn) denotes the probability that a raid 

from the base-aircraft location finds the i fcl1 target 

given that n aircraft survive enroute. 

( I>3 • 5 ) denotes the probability that any 

aircraft in a raid from the base-aircraft location 

survives until the k tJl release against the i*'* 1 target 

given that n aircraft are alive commencing the first pass, 

s t 

k = 1, or given that n aircraft are alive at the (K - 1) 
release for the second and subsequent passes, K = Z\ 

k = 1 ,2 ,3 f . . . f j . 

(D3.6) denotes the random variable which 

represents the number of the base-aircraft type alive 
at release on the k*'* 1 pass against the i^ target given 
that are assigned. The range of possible values 1st 
Xij k (Nij) = 0,1 ,2 , . . . ,N 1 j . For simplicity x 1 j k (N 1 j) is 
sometimes written as Xj^. 

III. DISTRIBUTION OF THE NUMBER OF AIRCRAFT 
ALIVE AT RELEASE ON THE FIRST PASS 

Under the above assumptions the event that x aircraft, 
x»-0, are alive at release on the first pass against a 
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target can occur In three mutually exclusive and 
exhaustive ways 2 



A raid of size N is undetected enroute, finds the 
target and x aircraft among N survive target defenses 
on the first pass, 

OR a raid of size N is detected enroute, but 
unengaged, finds the target and x aircraft among N 
survive target defenses on the first pass, 

OR a raid of size N is detected enroute, engaged, 
some n aircraft survive enroute, n == x, those n find the 
target and then x among n survive target defenses on the 
first pass. 

The extension of this type of reasoning used in 
conjunction with the laws of elementary probability 
theory produced the probability distribution of the 
number of aircraft alive at release on the first pass, 
Xiji, presented in Table I. 

IV. DISTRIBUTION OF THE NUMBER OF AIRCRAFT 
ALIVE AT RELEASE ON SUBSEQUENT PASSES 



The probability distribution of k = 2, was 

obtained from the distribution of by making the 

following observations l 

The stochastic process, £ X ijk J k " 

is §. discrete parameter finite Markov chain. That is, 

P [ X ijk ~ A k I "in " A 1 



I ^4 4 1 5 — 9 ° 



5 X 1 jk-1 “ x k-l 



■ p [ x i 



Jk z k I X ijk-1 ~ x k' 
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J 



The number of the J 
th 



base-aircraft type alive at 
release on the k 1 ' 11 pass against the 1 target Is less 
than or at most equal to the number alive on the (k-1) 
pass. That Is, 



st 



1 Jk 



“ X ijk-1 k " 



Therefore the conditional distribution of k - 2, 

given that ^=m, is binomial with parameters 

PHijk( m ) a ^d m as presented in Table II. 



By the theorem of total probabilities; 
N ^-J p lx 

ljk 



P Kjk “ *] " m fJ P [ X ljk - 1 | X 1 Jk-1 -■]*[* 



ijk-l 



= m] 



The synthesis of these facts yielded the probability 
distribution of X^j^* k - 2, presented in Table III. 



V. DISTRIBUTION OF THE TOTAL NUMBER 
OF PASSES AGAINST A TARGET 



The derivation of the probability distribution of 
necessitated one further assumption; 

(A3. 5) The number of a given base-aircraft type Alive 
at release on any pass over a given target is statistically 
independent of any other base-aircraft types alive over 
the same target on any pass. That is, X^j k is assumed 
independent of X^ y k / provided J f 

Under all foregoing assumptions the probability 
distribution of as presented in Table IV is derived 
in Appendix A. 
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PROBABILITY DISTRIBUTION OF X 
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elsewhere 
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elsewhere 
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CHAPTER IV 



THE PARAMETERS 
I. TARGET VALUE 

The value of a target may be characterized by a time 
dependent index*, i„e, 5 target value will change as the 
tactical situation changes over time* Furthermore, the 
relative values of targets in given tactical situations 
will depend strongly upon the doctrine adopted by the 
local military commander* If an index of target value 
is to provide a basis for decision, it must be formulated 
with consideration given to the whole social, political, 
economic, and ethical processes within which military 
action takes place. The intent of this section is not 
to formulate a scheme for target value indexing but to 
present some aspects of the problem which should be 
considered In any indexing plan 0 

Kaplan j2.j delimits areas of empirical study by 
which each element in the concept of military worth can 
be measured or appraised. These areas Include values, 
objectives, welfare, achievement, and worth, 

Kaysen £ 3 ] outlines criteria for selecting target 
systems asj 
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1. Military importance -- a rough classification 
of the value to enemy military operations of all types 
of equipment and supplies used by enemy forces J 

2. Percent direct (or indirect) military use -- the 
share of total output of a product or service which goes 
into military usej 

3. Depth — an Indication of the time available to 
the enemy for the organization of substitute consumption, 
alternate production, etc., before suffering military 
damage 

4. Economic vulnerability -- includes 

a. Ratio, of capacity to output, 

b. Substitutability for processes and equipment, 

c. Substitutability for product, 

d. Process and plant layout vulnerability, 

e. Recuperability J 

5. Physical vulnerability \ 

6. Location and size of the target system. 

In addition to the above criteria for indexing target 
value the effect on enemy morale should be considered. 

Hesse and Mitchell in * * Limited War Campaigns! 

Dacca Method* * [lo] have derived an index of target value 
which incorporates several of the criteria mentioned 
above. They consider military targets from a balanced 
ton standpoint using weighting factors for consumable and 
reuseable equipment, supplies, and manpower. This index 
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is quite involved^ but it provides a method of assigning 
target values for most conceivable military targets. It 
should serve well as a useful working tool for a first 
look at target value in the allocation problem. 

II. FUEL AND ORDNANCE LOADINGS 

To maximize damage to the enemy requires that air- 
craft be put over their targets with optimal ordnance 
loadings. A basic prerequisite to attainment of this 
objective is determination of fuel required to carry out 
the given mission. 

The following basic mission profile has been postulated 
to determine fuel requirements l 

1. Warm-up and take-off s 

2. Climb to cruise altitude and cruise to descent 
point for run-in to the target 

3. Descent to target run-in altitude 9 

4. Run-in to the target 9 

5* Ordnance deli very, 

6. Run-out from the targets 

7. Climb to cruise altitude and return to base 9 

8. Descent to the base for landings 

9® Landing and reserve. 

For ease of computation the following assumptions 
have been made® 

(A4.1) Aircraft carry full internal fuel and use 
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fully loaded 300 gallon external fuel tanks as required; 

(A4.2) A fuel weight of six and one-half pounds 
per gallon was used; 

(A4.3) No fuel Is burned and no distance over the 
ground Is covered in descents; 

(A4.4) All - external store racks for ordnance are 
equivalent on any particular aircraft type; 

(A4-5) There are no asymmetric load restrictions 
for launching aircraft; 

(A4-6) Five minutes fuel at normal rated power was 
allotted for warm-up and take-off; 

(A4-7) Thirty minutes fuel at sea level maximum 
endurance was allotted for landing and reserve. 

Utilizing the mission profile and the assumptions 
stated above one may formulate the basic fuel required 
equation as 

(4.1) Total fuel required = Take-off fuel + (2 x 

cruise fuel 'rate x 
(Range - Climb distance - 
Run-in distance)) + (2 x 
Run-in distance x Run-in 
fuel rate) + (Time over 
target x Target fuel rate) 

+ Landing and reserve fuel. 

After determining the total fuel required the next 
step is to ascertain if the fuel required is less than 
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the fuel available. Successful mission performance is 
constrained by 

(4.2) Fuel available - Fuel required H). 

Substituting for Fuel available in (4.2) one obtains 

(4.3) (Internal fuel + External fuel) - Fuel 
required ^0. 

Solving for External fuel results in 

(4.4) External fuel » Fuel required - Internal fuel. 

Once (4.4) has been solved it is possible to determine 

the number and size of external fuel tanks required and 
subsequently to determine the number of external store 
stations available for ordnance loading. 

The maximum ordnance load may be obtained by solving 

(4.5) Ordnance load = Maximum take-off weight - 

(Basic aircraft weight + 

Internal fuel + External fuel). 

Appendix B contains a program (PROGRAM 0RDL0AD) written 
in FORTRAN IV for the CDC-1604 digital computer which 
computes for a given aircraft mission the fuel required, 
the number of store stations available for ordnance, and 
the maximum ordnance loading. This program uses the 
assumptions and the algorithms outlined above for 
computations. Inputs required to the program, definitions, 
a program listing, and a sample output are given in 
Appendix B. 
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III. OTHER PARAMETERS 



The probability that an aircraft survives to make 

passes against a target is a function of the probability 

of detection by the enemy, the probability of engagement 

by the enemy, the probability of being killed enroute to 

the target by the engaging force , the probability of 

acquiring the target, and the probability of being alive 

thi 

at the weapons release point on the k pass against the 
target* Values for all of the above are required to 
obtain the distribution of W^. To minimize the objective 
function one also needs values for the probability that 

i.U 

the i target is killed by one pass given that b 
weapons are released on that pass. These probabilities 
will be considered as input parameters. 

Derivation of each of these probability distributions 
is beyond the scope of this paper 9 however , the remainder 
of this section will be devoted to discussion and/or 
references to derivations of these distributions. 
Throughout the investigation of these parameters , past 
data will in some instances be adequate and reliable, and 
the parameters could be estimated by Bayesian techniques. 
World War II , Korea, and Viet Nam experiences should 
result in extensive information which may be used to 
increase confidence in the estimates used. 



Koopman derives an expression for detection 
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probability both for visual and for radar detection. 
Electronic and other intelligence data coupled with 
characteristics of friendly radar should allow an 
estimate to be obtained for detection probability. Crout, 
Fay, and Harvey £ 5 ] have formulated a model describing 
a bomber’s penetration into hostile territory in the 
face of a given area defense strategy. Specifically an 
expression, is developed to compute the probability that 
an. attacking aircraft is killed before penetrating to a 
given depth. This model then accounts for the probability 
of detection, engagement, and kill. 

Simultaneously with passing through enemy area 
defenses friendly aircraft must attempt to acquire their 
assigned targets. Acquisition includes the processes 
of search, detection, identification, and flying the 
aircraft into a position to make a weapons pass against 
the target. Visual acquisition of the target is, of 
course, the most reliable and informative method of 
locating a given target. Due to circumstances such as 
weather, terrain, camouflage, aircraft speed, and aircraft 
altitude a pilot will sometimes have to rely on radar, 
infrared, microwave, or vectoring from personnel not in 
the aircraft to acquire his target. A combination of 
these methods should increase the probability an aircraft 
acquires its target. Erickson ^9j presents a study on 
the visual capabilities of a pilot searching for ground 
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targets,, Kamrass and Heckroth |l5j include a section 
on target acquisition which discusses detection 
probabilities and an analysis of types of sensors and 
methods of operation which will be most effective for 
reconnaissance and strike/ .missions-* 

Once an aircraft has acquired a target and is positioned 
for a weapons release pass it will usually be in the 
area of effectiveness of the point defenses of the target. 
These point defenses can run the gamut from no defense 
to a completely automatic system of surface-to-air 
missiles with a ground level to outer space capability. 

The probability an aircraft can successfully release 
weapons over the target is highly dependent upon the type 
of defensive weapon encountered. Selection of Aircraft 
for Tactical Air Missions ^24j discusses and gives values 
for probabilities of kill and the effect of various 
penetration aids associated with several types of point 
defense weapons. 

If an aircraft survives the point defenses 9 it will 
release weapons against the target. The probability of 
killing the target on one pass will be a function of the 
type and number of weapons released on that pass 9 the 
skill of the Individual pilot,, the weather,, the type of 
terrain,, and the type of target. The Naval Ordnance Test 
Station (NOTS) has worked on the development of weapons 
and the formulation of models which allow the probabilities 
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of kill of those weapons to be computed. NOTS technical 
publications include those of a probabilistic nature — 
Kusterer J ^8 , W.B, Simecka [l&] ? Verry [19] * K.D, 

Slmecka [26 , and Strang [28] J and those of a 

computational nature — Weldon and Young [l*+] * w » B. 
Simecka, et. al. [23 j , and * * Conventional Air-Delivered 
Strike Weapons’* [2?J . 

An alternative to derivation of the various 
probability distributions for all the parameters would 
be to use three estimates of each parameter, l.e., a best, 
an optimistic, and a pessimistic estimate. These will 
be referred to as BOP estimates. After obtaining 
solutions to the allocation problem using these BOP 
estimates, a sensitivity analysis should be performed 
to discover which parameters significantly affect the 
results* Following this analysis high confidence 
estimates for the most sensitive * ’parameters * * could 
be obtained using the distribution theory. 



28 



CHAPTER V 



CONCLUSIONS AND RECOMMENDATIONS 
I. CONCLUSIONS 

This analysis has demonstrated that the allocation 
of aircraft to fixed location targets is a complex 
problem,, If it had been reasonable to assume that 
attack aircraft make a single pass 9 as in strategic 
nuclear warfare 9 then the nonlinear model reported by 
Fiacco and McCormick j^22| could have been used 0 However 
in limited warfare 9 for which this model was developed, 
it is not reasonable to assume a single pass 0 Furthermore, 
we believe that a deterministic treatment of repeated 
passes would not suitably reflect the expected outcome 
in even an approximate fashion „ The substantiation of 
this assertion is dependent upon the generation of the 
probability distribution of 

The objective function for this model is non- 
dif ferentiable 9 this property restricts the class of 
usable and currently available nonlinear programming 
algorithms to that of Rosen £ 7 ] » 



29 



II. RECOMMENDATIONS 



It Is recommended that I 

1. This model be coded for use on an electronic 
digital computer utilizing the algorithm of Rosen |?J • 

2. The list of targets be grouped into type or 

character classes and an investigation made of the feas- 
ibility of using common values vfithin each group for the 
parameters PA^j(n), ) an d respectively. 

3. The list of targets be grouped into geographical 
classes and an investigation made of the feasibility of 
using common values within each group for the parameters 
PDij(n) and PE^fN,!!) respectively. 

4. A sensitivity analysis using BOP (best, 
optimistic, pessimistic) estimates of the parameters be 
made to determine those requiring further analysis. 

5. The feasibility of invoking the Central Limit 

Theorem on the expression W. = £ E X. be determined. 

1 j k 1JK 



30 



BIBLIOGRAPHY 



BIBLIOGRAPHY 



1. Operations Evaluation Group, Office of the Chief of 

Naval Operations, Department of the Navy. Search 
and Screening , by B, 0 * Koopman. 1946. OEG report 
no. 56 . 

2. The Rand Corporation. The Concept of Military Worth . 

by A. Kaplan. 7 May 19^. RM-37. 

3. The Rand Corporation. Notes on Some Historical 

As pects of Target Selection , by C. Kaysen. 

15 July 19 ^ 9 . RM-I 89 . 

4. The Rand Corporation. Optimal Tactics in Air 

Campaigns With Random Outcome's . by M. Dresher. 

February 19567 RM-1645, CONFIDENTIAL. 

« 

5. Naval Ordnance Test Station. A Probability Model 

for Bomber Penetratio n, by P. D. Crout , E. A. Fay, 
and J. R. Harvey. 12 August 1959* NOTS TP 2264. 

6 . Parzen, E. Modern Probability Theory and Its 

Applications . New York; John Wiley and Sons, Inc., 
I960;: 

7. Rosen, J, B. * ’The Gradient Projection Method for 

Nonlinear Programming. Part I. Linear Constraints * 
Journal of the Society of Industrial Applied Mathe- 
matics . v. 8, no. 1, March i 960 : 181-217. 

8 . Naval Ordnance Test Station. The Expected Number of 

Weapon Deliveries Needed to Damage a Target . by 
D. F, Kusterer. 3 April 1961. NOTS TP 2658, 
CONFIDENTIAL. 

9. Naval Ordnance Test Station. Air -to -Ground Visual 

Acquisition of Tactical Targets . by R. A. Erickson. 

3 November 1961 . NAVWEPS report 7804, NOTS TP 2801, 
SECRET. 

10. Laboratores for Applied Sciences, University of 

Chicago. Limited War Operations . December 1961 . 
LAS-TR-306-29, SECRET. 



31 



11. Roser^ J«B. 9 9 The Gradient Projection Method for 

Nonlinear Programming. Part II. Nonlinear 
Constraints. 9 9 Journal of the Soci ety of, Industrial 
Applied Mathematics . v. 9s> no. 57~December 19^1. 
51^-532. 

12. Parzen, E. Stochastic Processes . San Francisco. 

Holden-Day j,“inc. , 19627" 

13 . The Rand Corporation. Recent Developments in Non- 

linear Programming , by P„ Wolfe"."** May 19^2 . R-4oi-PR. 

14. Naval Ordnance Test Station. A Handbook on the 

Blast Effectiveness of Conventional Warheads Against 
Limited-War Ground Targets . by R.G. Weldon and G. H. 
Young. June 19^2. NAVWEPS report 7912 s NOTS TP 
2923 9 CONFIDENTIAL. 

15. Cornell Aeronautical Laboratory, Inc. Evaluation 

of Tactical Air Weapons Systems . Volume II s> by M. 
Kamrass and E. E. Heckroth. September 19^2. CAL 
report GM-l496-G-5s> SECRET— RESTRICTED DATA. 

16. Naval Ordnance Test Station. Delivery Accuracy of 

Air-Launched Weapons Delivered Against Ground Targets 
In An Unfamiliar Environment , by W„ B. Simecka. 
December 19&>. NOTS TP 3148;, CONFIDENTIAL. 

17. Graves j, R. L. and Wolfe, P. Regent Advances In 

Mathematical Programming . New YorkS McGraw-Hill 
Book Company j, Inc. 9 19o3 » 

18. Dorn, W. S. 5 '’Nonlinear Programming =~ A Survey 99 

Management Science , v. 9, no. 2, January 1963* 

171=208. 

19* Naval Ordnance Test Station. Round Requirements : t to 

Damage a Target At A Specified Level , by W.R. Verry. 
May 1963." NAVWEPS report 8362, NOTS TP 3258. 

20. Research Analysis Corporation. Programming Under 

Nonlinear Constraints by Unconstrained Minimi zatlon S 
A Primal -Dual Method , by A. V. Fiacco and G. P. 
McCormick. September 1963. RAC TP 96. 

21. Cornell Aeronautical Laboratory, Inc. Tdctlcal Air 

Warfare Force Allocation . October 1963. SECRET. 

22. Research Analysis Corporation. Algorithm for Non- 

linear Programming " Re source -Allocation Methodology , 
by A. V. Fiacco and G. P. McCormick. November 1963. 
RAC TP 108. 



32 



23 



. Naval Ordnance Test Station, Combat Effectiveness 

of Alr-to-Surfaoe Weapons with Conventional Warheads, 
by W, B, Slmecka and others. November 1963, NAVWEPS 
Report 8393, NOTS TP 3308, SECRET. 

24. Arthur D. Little, Inc. Selection of Aircraft for 

Tactical Air Missions . 15 November 1963 . SECRET. 

25 . Hadley, G. Nonlinear and Dynaml c Programming . 

Palo AltoS Addison-Wesley Publishing Company, 1964. 

26. Naval Ordnance Test Station, Fractional Single-Shot 

Kill for A Cluster Weapon Delivering A Circular 
Pattern Against A Rectangular Target Area , by K. D. 
Slmecka. January 1964. NAVWEPS report 8402, 

NOTS TP 3323, CONFIDENTIAL, 

27. Naval Ordnance Test Station, Conventional Air- 

Dell vered Strike Weapons . 1 January 1964. NOTS 
TP 3408, CONFIDENTIAL* 

28. Naval Ordnance Test Station. Warhead Lethality 

Analysis and Program Computations . by J. Strang. 
February I 964 . NAVWEPS report 8468, NOTS TP 3424, 
SECRET. 

29. Center for Naval Analyses. Planning Factors for the 

Analysis of Tactical Air Force Levels . by G. Haering. 
May 1964. CNA study S 77 , SECRET. 

30. Research Analysis Corporation. The Sequential Un- 

constrained Minimization Technique for Nonlinear 
Programming . Algorithm II . Optimum Gradients by 
Fibonacci Search , by A. V. Flacco and G. P. McCormick. 
June 1964. RAC TP 123. 



33 



APPENDIXES 



APPENDIX A 



DERIVATION OF THE PROBABILITY 



DISTRIBUTION OF W^I^) 



The event = 0 occurs if and only if ^ *ijk = or 
if and only if X^.^ = 0 for all j and all k. It was noted 
that X^j k — X^j k _^, k = 2,3*4, ...,R^j, 30 that X^j k = 0 
for all j and all k if and only if X^^ = 0 for all J . 

It was assumed that X^j k is independent of X^/ k / provided 
that j ^ j 7 , so thatS 

P [w t = o] = p[E Z X 1Jk =0]=p[x iJk = 0 for all j and all k] 

J k 



= p[x = 0 for all 



jl = n p[) 
J J 



ijl 



= o] 



The distribution of 1 was obtained in a 

slightly different manner. The event = 1 occurs if 
and only if for exactly one j, X^j^ = 1 and X^^ = 0> 
but the event = 2 occurs if for exactly one j , 

= 2 and X^j 2 = °* or for exactly one j, X^^ = l»X 1 j 2 
and X^j^ = 0, or if for exactly two j*s, X^j^ = 1 and 
X^j 2 = 0* This enumeration can be continued and was up 
to =4, it produced a lengthy table from which a 
recursion relation was developed. The use of this relation 
in conjunction with the probability distribution of X^ k 
produced the probability distribution of in the following 



manner 
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(D6.1) Let faco denote, for all combinations of. Where 
* ’combinations * 9 is used in the sense of combinatorial 



analysis, e.g., Parzen f^6 

p h - 1] - *[p x in = i 



, then 



P [ X 1J2 = 0 I X ijl - P [ X 



ini 



P [ W 1 - 2] - = p [*i» - 2-J p[x ij2 - «! K 



J a x =0 



J1 ~ 2 " a l 



• P [ X 1J3 I X 1J2 = “J n * 3 P [ X 1„1 - °] 

+ ftco p = X ] p [ X U,2 = ° 1 X UU = l J 

J i > J 3 

* P [ X lJal =1 ^ P [ X ij »2 = 0 l X ij a l = X ] , P [ X inl = °] 



note: p [x 1Jk = x 1 x 



ijk-l 



= m] = 0 



for all x > m and 



p [ x lJk ’ 0 I X ljk-1 ' °] * Table «• 

P ^ = 3 J “ f ajo S P l^n “ 3 ' 



,P [ X iJ2 a l " a 2 I X ijl 3 “ a J P [ X ij3 “ ”2 1 "iJ2 " "1 



= a, X, ., = a. - a 



' * P [ X i J4 0 I X iJ3 " a 2-l P [ X inl " °] 



(continued on the following page' 
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+ 



+ 



2 3-a,-l a,-l r -| 

E E 1 X E P I X. . , = 3 - a - b J 

=1 b i:L .o b 21 =o L lJl1 1 11 



E 

faco a, 

Ji»Ja 



P [ X iJx2 = b ll 1 X iJxl 



-3 - a l - 



X iJx3 0 I X Ui2 " b lJ 



P [ X iJal “ a l “ b 2l] P [ 



X ija2 = b 21 I X ij a l 



= a l " b 2ll 



P [ X iJ.3 = - ? 1 X «,2 = b ^n2j,,0. 



p [ x lnl “ °] 



faoo P = ^ = ° I 



'A 



Jx , Ja » J 3 



P [ X i Jal - '} P Kj.2 - °.| X lJ a l ■ *] 

p [w^] P [ X Ua2 “ ° I X 13al “ p [ x inl “ 0 
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By Induction the general term of the probability 
distribution of was obtained and is presented in Table 
IV, where all sums are defined in the positive sense, 
e 

i.e., the sum E f(y) is defined if and only if e - d 
y=d 

and the following definitions were usedS 
(D6.2) R denotes the max (R, ,) 

i.J J 

(D6.3) r denotes the max (R S J) 



Summation Operators l 



(D6.4) 


r=l 

B, = TT 

1 k i =2 


b ik. -I” 1 
X E 

b i ki = 0 




(D6.5) 


J-l 

C . = TT 

J k=2 


a k =J.k 




(D6.6) 


W— 1 

A, - E 

1 a^O 


r : x “k-i; 1 

k=2 a k=0 




(D6.7) 


w-1 

k 7 = E 

af 1 


w-a, -1 a, -1 

’x B. H 

b ll = 0 d 21“ C 


1 B * 


(D6.8) 


w-1 

A, = E 

a l= j-1 




J-l a. n -a t -l 1 

tt t=i E b B. 

4*1 b tl =0 


B, 

b jx =0 J 




j ~ 3 p • • • 


j,J and a Q - w 
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Summation Operands I 

(D6.9) O tm (l) = P [x 1Jtl = a t . 1 - a t (l-4 tm ) - b tl ] 
• P [ X ij t 2 = b., - b 



tl “M | X lJ t l a t-l " a t* 1_d tm* " b tl] 



r-2 



V 2 



" r 

TT P [ 



= b. , - b 



U t k t +1 tk^. tk^+l | ij^k*. tk^-1 u tk 



= b. 



- b. 



t“t 



] 



* P [ X iJ t r b tr-l I X iJ t r-l ~ b tr-2 " b tr 

• P [ X iJ t r + l = 0 I X iJ t r = b tr 



J 



-J 



where d^ = < 0 if i ^ J 

1 if i = j 



r-2 



(D6.10) P. n = Z tt P f 
1X Ji k=l L 



X ij x k " a k- 



-1 ” a kl 

J 1 

[ X ij A k+l = a k " a k+l I X ijxk = a k-l " a k] 

[ X lUr = a r-l I X i^r-1 = a r-2 ~ a r -J 
iJir+1 



P [x tl = 0 



X ijir - a r 



-J " P [ X 



n^j; 



ini 



-o] 



(D6.ll) P la = S 



m 



faco 



n G tm (1) 
t=l tm 



tT 



Jl I • •• t 



+ y - - ~ j 

P [ X lnl ■ °] 



V 

m = 2 ;3 i ; ; ; *J-1 and J ^ J: 



(D6:i2) T ±j = n G tJ (i) J - 2 

t = l 
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Therefore under the above definitions, given values 



for R and J, the probability distribution of 
in Table IV can be expressed in closed form for 
computational purposes. 



presented 
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APPENDIX B 



PROGRAM ORDLOAD 

The FORTRAN IV program to compute the fuel and 
ordnance loadings is presented in this Appendix. A 
brief description of the imput data, the definition of 
variables, a program listing, and a sample output will 
be given. Data for the sample problem was selected at 
random, and the resultant output bears little semblance 
to any realistic situation that might occur. The program 
has been dimensioned to handle six aircraft type, six 
bases, and fifty targets. These dimensions may be 
expanded as necessary to suit a particular user’s need. 
Common blocks have also been established to facilitate 
incorporation of this program as a sub-routine in a 
larger allocation problem. 

The data input is presented below. Variables with 
subscript (I) require one data cqrd for each aircraft 
type. 

Card i: FORMAT (3110, F10. 0 ) l 

Fields l-io: NTYPE = the number of different 

aircraft type, 1 £ NTYPE 6j 

Fields 11-201 JBASE = the number of friendly 
bases, 1 - JBASE - 6J 
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Fields 21-30 l KTGT = the number of targets s 
1 - KTGT ~ 50 ; 

Fields 31-40: TIME p* the time ( tenths of 

hours ) over the target for delivery of ordnance. 

Card 2,...,NTYFE+i: FORMAT ( 5FX0. 0,110 )l 

Fields 1-10: GAXTOW(I) = maximum gross aircraft 



take-off weight (pounds )£ 

Fields 11 -20: BACFTW(X) = basic aircraft weight 

( pounds ) J 

Fields 21-30: ENTFW(X) = maximum internal fuel 



capacity (pounds): 

Fields 31-40: EXTFW(X) = maximum external fuel 

capacity (pounds): 

Fields 4X~50: GAXBMW(I) = maximum gross ordnance 

load (pounds): 

Fields 51-6 o: XRACK(I) = maximum number of 

external store racks on an aircraft: 



Card NTYPE+2,...,(2xNTYFE)+i: FORMAT (5F10.0): 

Fields 1-10 : FTAKE(X) - fuel (pounds) required 

for warm-up and take-off: 

Fields ll-20: FCLIMB (X ) = fuel (pounds) 



required for climb to cruise altitude: 

Fields 21-30: FLOXT(I) = loiter fuel rate 

(pounds per minute)-- not required in the current 
formulation of the program: 
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Fields 31-40! DASHSP(I) = speed (knots) over 
the target at military power; 

Fields 41-50: FLAND(I) = landing and reserve 

fuel (pounds). 

Card (2xNTYPE)+2,. , (3xNTYPE)+i: FORMAT (2F10.0): 

Fields 1-10: DISTML(I) = distance (nautical 

miles) covered in the climb to cruising altitude ! 

Fields 11-20! DISTRI(I) = run-in distance 
(nautical miles) to the target. 

Card (3NTYPE)+2! FORMAT ( 2F10.0)! 

Fields 1-10! ALT1 = cruise altitude (thousands 

of feet); 

Fields 11-20! ALT2 = run-in altitude (thousands 

of feet); 

Card (3xNTYPE)+3,..., (4xNTYPE)+2! FORMAT (4F10.0)! 

Fields 1-10! SFC(I)„= sea level cruise 
specific fuel consumption (pounds per nautical mile = 
ppnm) ; 

Fields 11-20! SFCMIL(I) = sea level military 
power specific fuel consumption (ppnm); 

Fields 21-30! SFCK(I) = cruise specific fuel 
consumption altitude correction factor (ppnm per thousand 
feet altitude); 

Fields 31-^0! SFCMILK(I) = military power 
specific fuel consumption altitude correction factor 
(ppnm per thousand feet altitude). 



42 



Card (4xNTYPE) + 3 » . . . , (4xNTYPE)+2+( JBAS£):cKTGT) • : FORMAT (FIO.O): 
Field 1 - 10 ! RANGE(J,K) = range (nautical miles) 
from Jth base to the Kth target. 

Definitions of other variables used in the program 

are! 

FCRUSA(I) = cruise altitude fuel consumption rate 

(ppnm) ; 

FCRUSL(I) = run-in altitude fuel consumption rate 

(ppnm) ; 

FDASHM(I) = run-in altitude military power fuel 
consumption rate (ppnm)! 

FAVAIL(I) = maximum total fuel available for 
aircraft type I! 

FRQRD(I , J,K) = fuel required for aircraft type I 
to complete mission from base J to target K and return! 

EXTFRQ(I, J,K) = external fuel required for air- 
craft type to complete mission from base J to target K 
and return! 

KRACK(I,J,K) = the number of external store racks 
available on aircraft type I to target K from base J! 

0RDLD(I 5 J„K) - ordnance load (pounds) aircraft 
I can carry to target K from base J. 

A complete program listing and a sample output, which 
uses the data at the end of the program listing, are given 
on the last four pages of this Appendix. The sample 
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output tables the aircraft number, the base number, 
the target number, the range from base J to target K, 
the total fuel required (lbs.), the Internal and external 
fuel required (lbs.), the number of ordnance racks available, 
and the ordnance load capability (pounds). 
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